Optimal purification of single qubits 
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We introduce a new decomposition of the multiqubit states 
of the form p® N and employ it to construct the optimal single 
qubit purification procedure. The same decomposition allows 
us to study optimal quantum cloning and state estimation of 
mixed states. 

PACS Nos. 03.67.-a, 03.65.Bz 



Quantum data processing [||, such as quantum cryp- 
tography, teleportation or computation, rely on an ample 
supply of qubits in nearly pure quantum states. However, 
maintaining unknown and uncorrupted entanglement or 
superpositions in the presence of the coupling to the envi- 
ronment is usually a difficult task and a great deal of re- 
cent research has been focused on issues such as entangle- 
ment purification J2|. In particular, several entanglement 
purification protocols have been devised but so far none 
of them has been proved to be optimal; that is, it is not 
clear what is the maximal fraction of purified entangled 
particles for a given initial state and a prescribed degree 
of purification. This paper addresses a less complicated 
question of purifying unknown quantum superpositions 
of qubits; however, it provides a definite and construc- 
tive answer regarding the optimality of the single-qubit 
purification procedures. The solution of the problem in- 
volves a particular decomposition of density operators of 
multiqubit systems. This decomposition seems to be a 
useful tool on its own as it also gives some insight into 
the related problems of the optimal cloning || and the 
optimal state estimation of mixed states. 

Suppose we are given N identical qubits, each in an 
unknown state described by a density operator 

p = ^(l + Afi.^)=ci|lfl)<l f t|+c |0 f t)(0 fI | > (1) 

where An is the Bloch vector, |n| = 1, c\ — cq = A > 0, 
ci + Co = 1 and |0#), are the two eigenstates of p. 
The goal is to obtain the maximum number of qubits in 
a state "as close as possible" to the unknown state 
i.e. to purify the state p. Before specifying more precisely 
the problem, let us re-write p as a convex sum 



P = A \ (1 + ft •<?) + (!- A) h 



(2) 



maximally mixed state 4l. This decomposition implies 
that given a set of N identically and independently pre- 
pared qubits one can distill from the set at most AiV of 
perfectly purified qubits. In the following we show that 
this limit is achievable when ./V — > oo and we describe 
the optimal purification procedure for any finite N. 

In general, a purification procedure V consists of a set 
of operations and measurements on the TV qubits and 
perhaps on additional ancillas. Depending on the out- 
comes of the measurements some of the qubits will be dis- 
carded and we will end up with M < N qubits whose re- 
duced density operators are closer to the ideal state 1 1^) . 
Thus V is described by a set of linear completely positive 
maps Vm such that Vm{,P® N ) = Pm where pu is an un- 
normalized density operator describing M qubits. The 
probability of obtaining M qubits is Pm = ^{pm) and 
their corresponding fidelities F M = (l,t\pk\^n), where 
Pk is the normalized reduced density operator of qubits 
k = 1,2,..., M. These probabilities and fidelities depend 
both on the direction n and the length A of the Bloch 
vector of the original qubits, and thus we write them as 
Pm(A, ft) and F M (X, n). Given that we are not interested 
in the fidelity of a particular qubit, nor in a particular 
direction of the Bloch vector but rather in their average 
values, we define 



(3a) 



P M (X) = / dnP M (\,ft), 



A I 



where the first term on the r.h.s. represents the puri- 
fied state | (1 + ft - a) = and the second one the 



Fm {x) = Pm\\) J M Pm (a ' h) Ji £ F ^ (A ' " } ' (3b) 

where the integral refers to the complete solid angle. 
Thus, any given V is characterized by these two quanti- 
ties namely, the probability of obtaining M valid qubits 
and their corresponding fidelity. This fact implies that 
there is no unique way of defining optimality. For ex- 
ample, a procedure might give with small probability a 
large number of qubits of high fidelity, whereas another 
one might give with a large probability a smaller num- 
ber of qubits of higher fidelity. Nevertheless, it turns 
out that there is a procedure V which has the follow- 
ing unique properties: (1) given any other purification 
procedure V' characterized by P' M and F' M , it is always 
possible to obtain these values by first applying V and 
then some prescribed operations and measurements; (2) 
after applying V, it is impossible to increase the fidelity of 
the outcoming qubits even at the cost of decreasing the 
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number of qubits. The first property is a consequence 
of the fact that V is, in a certain sense, reversible. It 
implies that after applying it, we can still obtain the op- 
timal results whatever we mean by optimal. The second 
property implies that after applying V, it is impossible 
to purify further. As a consequence, our procedure is op- 
timal for any reasonable definition of optimality. On the 
other hand, we introduce now an important concept that 
will be used later on to prove optimality of our proto- 
col. We say that a procedure V s is symmetric if: (i) the 
reduced density operators of the M qubits are identical; 
(ii) the maps V S M are covariant , that is, 



U® M V a M (p 6 



(C/t) 



(4) 



for all U € SU(2). These conditions imply that the prob- 
ability of obtaining M qubits is independent of n and that 
the fidelities of all the qubits are the same, which sim- 
plifies the calculations. Moreover, they are the only ones 
that have to be analyzed since for any arbitrary purifica- 
tion protocol V one can construct a symmetric one which 
gives exactly the same values of Pm(A) and Fm(A). In 
order to do that, we proceed as follows: (1) we apply the 
same random unitary U operator to the N qubits; (2) 
we apply the procedure V; (3) we apply the operator [/' 
to the M resulting qubits; (4) we permute the M qubits 
randomly. 

We shall start with a set of N = 2J qubits (the 
2 J + 1 case can be solved in an analogous way) and de- 
compose their state p® N into the sum of density oper- 
ators pj a which are supported on orthogonal subspaces 
Sj, aj = span{|j, m, aj); m = -j, -j + 1, . . . , j}, where 
j — 0, 1, . . . , J, a>j — 1,2,..., dj, and 



1,2, 



2J 
J-j-l 



C? * J), 



(5) 



dj = 1. States \ j, to, a) are constructed from an arbitrary 
basis {|0), |1)} of a single qubit in such a way that: (i) 
for a fixed j and a, vectors \j, to, a) form a basis for the 
(2j + 1)— dimensional irreducible representation of SU(2); 
(ii) for a fixed j and to, vectors \j, to, a) form a basis 
for the dj— dimensional representation of Sn (symmetric 
group) corresponding to the Young diagram (2j, 2J— 2j). 
In particular, we choose 



|j,to,1) = \j,m) ® |*-) 



iJ-j 



(G) 



where |j, to) is the symmetric state of 2j qubits with j— to 
qubits in the state |0) and j + to qubits in the state |1), 
and |^_) is the singlet state of two qubits. The rest of the 
states are constructed by using permutation operators, 
IT £ Sn, 



\j, TO, a) =^ h i U i\j> m > !) 



(7) 



and the coefficients hi are chosen so that the correspond- 
ing states are orthonormal. Although we have defined 



the basis \j,m,a) in terms of the states \i) (i = 0,1), 
we could have used any other basis — Ufi\i) where 
Un 6 SU{2) to define \j,m,a)a. One can easily check 
that 



\j,m,a)ji = U? N \j,m,a) 



E^,mWli.*»'.«). (8) 



where D 3 m , m (n) are the corresponding matrix elements. 
Using (||) one immediately sees that the subspaces Sj t0l 
are the same if we take any basis \j, to, a)^ to construct 
them. We also define a set of unitary operators Uj. a 
which transform Sj Q into Sj,i while leaving the other 
subspaces untouched. We take Uj jCt \j,m,a) — \j,m,l) 
and Uj a \j, to, 1) = \j,m,a). Again, using (|J) one can 
easily check that these definitions are valid for all n, i.e. 
Uj, a \j,m, l)n = \j,m,a)a. 



The decomposition itself reads 



where 



j=0 J a—1 



2] + l _ 2j+l 
c 



Pi = dj[c ci]- 

Cl-CQ 

Pj-a = U 3,aPjsU}, a for Q ^ 1, 
PJ,1= Pi ®(l*-)(*-|)® J - J . 

i 



(9) 



(10a) 

(10b) 
(10c) 



Pi = 2j+i 2j+i E c o m(J i +m \h m )n{j,m\- (lOd) 



This last density operator can be conveniently written as 
convex sum of density operators corresponding to iden- 
tical pure states of 2j qubits: 



Pi = 2 J + 1 _ 2j + l ( 2 J + !) 



dn 



(ii) 



where n(9) = c\ cos 2 (6»/2) + c sin 2 (6>/2) and 

|^(M)}=^^#|l}, + ^^#e^|0), (12) 
^n(6») V n (^) 

Any state of the form p® N can be decomposed into the 
Pj t a components as in Eq. (j^); this can be seen by writ- 
ing p® N = J3/s=o c o c i~ k Pk, where Pk is a projector on 
the subspace with exactly k qubits in state |0)^. Ex- 
pressing this projectors in terms of the basis |j, to, a)fj 
immediately gives (|9|). 

Let us now use decomposition (^) to introduce our pu- 
rification procedure. Given the state p® N we perform a 
measurement defined by the projections on the mutually 
orthogonal subspaces Sj t0l . This provides us with one of 
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the states pj. a . If a ^ 1 we apply the appropriate unitary 
transformation uj a pj, a Uj,a to obtain p^\. The result- 
ing density operator, — pj ® (l^-) |)® J ~-' gives 
a clear separation between the first 2j qubits in state 
Pj and the remaining qubits paired into J — j singlets. 
Since the singlet state carries no information about state 
p, we can discard the last 2{J — j) qubits. This procedure 
gives, with probability Pj, 2j qubits in state pj. The re- 
duced density operators of the resulting qubits coincide 
and have a fidelity 



1 



(2j + l)c 



2J + 1 



Cl 



„2i+i 



2i+l 



ci - c 



(13) 



Thus, we have that the purification procedure gives [cf. 
Eq. (||)] P M = p j=M /2 and F M = fj=M/2 if M is even 
and Pm = otherwise. 

In order to illustrate our considerations with an ex- 
ample, let us take the simplest case of N = 2 qubits. 
The density operator p® 2 can be decomposed into two 
terms supported respectively on the singlet (antisymmet- 
ric, /So,i) and the triplet (symmetric, Si i) subspace, 



p® 2 =J2PiPj,i=coCi\*-)(*-\ + 



(14) 



3=0 



+ ( C 2|0)(0| + Co c 1 |* + )(* + | + c?|l)(l|). 

Our purification procedure consists of projecting onto the 
symmetric or antisymmetric subspaces. After success- 
ful projection on the symmetric subspace, we keep both 
qubits; after successful projection onto the antisymmet- 
ric subspace, we discard both qubits since the singlet 
state does not depend on ft, and therefore the qubits do 
not carry any information about the original state. This 
procedure is symmetric and gives ^(A) = 1 — coCi = 
(3 + A 2 )/4 = 1 - P (A) and F 2 (A) = ci(l - c /2)/(l - 
cqCi) > ci. Thus, with probability P 2 (A) we obtain a 
higher fidelity than the original one. 

In general, the average distillability factor (or yield) 
and the mean fidelity of our procedure are given respec- 
tively by 



d n (\) = YtViTn - A + + (){ L/v " ) - ,ir,a) 



F N (\) = j^Pjfj ^ 1 " + ° (1/Ar2) - ' ' r>h] 



3=0 



In the limit N ^> 1, the mean fidelity tends to one 
whereas the yield tends to A, the length of the Bloch 
vector. In fact, one can check that in this limit the dis- 
tribution pj becomes narrower and narrower (that is, the 
width divided by the mean value tends to zero). This 
means that if N is sufficiently large we will basically ob- 
tain iVA qubits in the unknown pure state \l)n- Thus, 



in the limit N — > oo our purification procedure gives the 
maximum yield. The question is: how good is this pu- 
rification for finite N? In fact it is the optimal one! 

In order to show that we have an optimal procedure 
we will first show that given any other purification pro- 
cedure V it is always possible to obtain the same val- 
ues P' M and F' M once we have applied our procedure V . 
Suppose we have applied V to the N qubits, obtained the 
outcome corresponding to a particular subspace Sj >a , ap- 
plied the corresponding operator Uj , and discarded the 
last 2( J — j) qubits. In order to obtain the results of the 
procedure V we just have to add J — j pairs of qubits in 
singlet states, apply the operation Uj, a to the whole set, 
and then apply the procedure V. It is straightforward 
to show that this whole procedure gives the same values 
of P' M and F' M . The intuitive idea behind all this is very 
simple: according to (||) the density operator p® N can 
always be regarded as describing the situation in which 
we have prepared with probability Pj/dj the state pj >a . 
Thus, our measurement only identifies the state which 
has been prepared without disturbing it, whereas the op- 
eration Uj t0l is unitary. Therefore, this procedure is to all 
practical purposes reversible. 

Now, we complete the proof of optimality by showing 
that after having applied V, it is impossible to increase 
the averaged fidelity of the outcoming qubits. The idea is 
to show that any purification procedure which produces 
1 qubit out of 2j qubits in state pj has at most a fidelity 



-fi(A) = fj given by (pL3|) . This automatically implies 
that there exists no procedure that produces M qubits 
(for any M) with fidelity higher than fj for if it existed, 
we could construct a procedure that gives 1 qubit of fi- 
delity larger than fj by simply discarding the remaining 
M — 1 qubits, which contradicts our previous statement. 
Therefore, we restrict ourselves to symmetric procedures 
that produce one qubit out of 2j qubits in a state pj . Let 
us consider a general linear and completely positive map 
V\. It follows from the covariance condition (^) that for 
any single qubit in state |$) we have 



Pi \mmT 23 \ =x\y)(*\+ y \^)(y- 



(16) 



where x,y > 0, x + y < 1 and (>]/ ± | , I') = (this can be 
easily proved by showing that the result must commute 
with any rotation along the axis defined by the Bloch 
vector corresponding to | V I')). Using Eq. ( |ll| ) and the 
fact that Pi is linear we can obtain Pi(pj) and express 
the fidelity of Pi(pj) /ti[Pi(pj)] in terms of a; and y. The 
maximum value occurs for (x,y) — (1,0) and gives ex- 
actly fj, which completes the proof. 

Interestingly enough the methods introduced above 
can be used to study the optimal cloning H[j]|| of mixed 
states. In this case the goal is to start with N qubits all in 
the state p (Q), and produce exactly M > N qubits each 
in a state as close as possible to the state 1 1) T -j M. More 
specifically, defining the average fidelity as in (|3b|), we 



3 



define the optimal cloning procedure C as the one which 
gives a maximal value of Fm- Again, one can easily show 
that symmetric procedures C s are optimal by using the 
same four steps as before. Moreover, since to all practical 
purposes the purification procedure defined above is re- 
versible, we can study independently the optimal cloning 
of 2j qubits in the state pj and at the end average with 
the corresponding probabilities pj . Let us then consider 
a symmetrical procedure that produces M qubits out of 
2j qubits and that is described by a linear trace pre- 
serving positive map Cj . The corresponding fidelity will 
be (lft\trM-i[Cj(pj)]\la) , where trM-i denotes the trace 
over any M—l particles. Using the property (^]) we have 



tr 



M-l 



where 



x — r 2j,M 



M(2j + l) + 2j 
M(2j + 2) 



(17) 



(18) 



since otherwise Cj would give for pure states better clones 
than the optimal cloning procedure Using the de- 
composition (p"l|), linearity, and (17) we obtain that the 
optimal value occurs for the largest value of x, i.e. for 
x = F^j M' Thus, the fidelity of the optimal cloning pro- 
cedure for mixed states is 



TTTIllX 

r N.M 



J 

3=0 



ppur 
r 2j,M 



Si + (i 



^■mXW,) , (19) 



where p 3 , f. 



and F 2^M 



are given m 



@ respectively. Defining XT*, = 2F mix 



Xn x m = ^%(M + 2)/M 



where 



|0a|), flljj), and 

— 1 we obtain 



N.AI 



acceptable input states for the optimal quantum don- 
ers. The decomposition, apart from having applications 
in restoring quantum superpositions via the single qubit 
purification, seems to be a useful mathematical tool on 
its own. 



\ mix 

A N,oo 




N 100 

FIG. 1. Maximum achievable length of the Bloch vector 
-^jvJoo in the cloning N — > oo as function of the initial number 
of copies of the state p. The values of the initial lengths of the 
Bloch vector are, from bottom to top, A = 0.2, 0.4, 0.6, 0.8, 1. 
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problem of optimal measurements. 
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(20) 



with -F 1 ™^, is the maximal fidelity in the N — > oo cloning 
of mixed states. Actually, using the same arguments as in 
H one can show that this fidelity is the maximal fidelity 
with which one can estimate the state |1^) starting from 
N copies of the state p. This way of estimating copies 
of mixed states might be interesting in the following con- 
text: Alice is sending N qubits in an unknown pure state 
to Bob via the depolarizing channel and Bob's goal is to 
determine the state that Alice is sending. In Fig. 1 we 
have plotted A™^, as a function of N for various values 
of the length of the Bloch vector. 

In conclusion - in this paper we have introduced a new 
decomposition of the multiqubit states of the form p® N . 
This decomposition was instrumental in constructing the 
optimal single qubit purification procedure (and proving 
that it is indeed optimal) and in extending the remit of 
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Note that this scenario is not equivalent to the problem 
of generating M output states as close as possible to the 
mixed state p. 
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